Edge of chaos of the classical kicked top map: Sensitivity to initial conditions 
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We focus on the frontier between the chaotic and regular regions for tlie classical version of 
the quantum kicked top. We show that the sensitivity to the initial conditions is numerically well 
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characterised by ^ = eg' , where = [1 + (1 — q) x] i-' (ej = e^), and \q is the g-generalization of 
the Lyapunov coefficient, a result that is consistent with nonextensive statistical mechanics, based on 
the entropy Sq = {1 — Yl^pD/iq ~ 1) (5'i = — Y2iPi Inpi)- Our analysis shows that q monotonically 
increases from zero to unity when the kicked-top perturbation parameter a increases from zero 
(unperturbed top) to ac, where Qc — 3.2. The entropic index q remains equal to unity for a > a^, 
parameter values for which the phase space is fully chaotic. 



I. INTRODUCTION 

Non-linearity is ubiquitously present in nature, e.g., fluid turbulence [1], extinction/survival of species in ecological 
systems [2], finance [3], the rings of Saturn [4], and others. Consistently, the study of low-dimensional non-linear 
maps plays a significant role for a better understanding of complex problems, like the ones just stated. In a classical 
context, a main characterisation of the dynamical state of a non-linear system consists in the analysis of its sensitivity 
to initial conditions. From this standpoint, the concept of chaos emerged as tantamount of strong sensitivity to 
initial conditions [5]. In other words, a system is said to be chaotic if the distance between two close initial points 
increases exponentially in time. The appropriate theoretical frame to study chaotic and regular behaviour of non-linear 
dynamical systems is, since long, well established. It is not so for the region in between, edge of chaos, which has 
only recently started to be appropriately characterised, by means of the so-called nonextensive statistical mechanical 
concepts [6]. In this article we study the sensitivity to initial conditions at this intermediate region for the classical 
kicked top map and its dependence on the perturbation parameter a. The sensitivity to initial conditions is defined 
through ^ {t) = lim||Ar(o)||-»o [[^^|o||| , where Ar (t) represents the difference, at time t, between two trajectories in 
phase space. When the system is in a chaotic state, ^ increases, as stated previously, through an exponential law, i.e., 

^(i)^e^^* (Ai>0), (1) 

where Ai is the maximum Lyapunov exponent (the underscript 1 will become transparent right a- head). Equation 
(1) can also be regarded as the solution of the differential equation, ^ = Ai ^ . In addition to this, the exponential 
form has a special relation with the Boltzmann-Gibbs entropy Sbg = ~ ^iPi In Pi- Indeed, the optimization of Sbg 
under a fixed mean value of some variable x yields p{x) oc e~^^ , where P is the corresponding Lagrange parameter. 

Excepting for some pioneering work [7], during many years the references to the sensitivity at the edge of chaos 
were restricted to mentioning that it corresponds to Ai = 0, with trajectories diverging as a power law of time [5]. 
With the emergence of nonextensive statistical mechanics [6], a new interest on that subject has appeared, and it 
has been possible to show details, first numerically [9] [10] [11] [12] [13] and then analytically (for one-dimensional 
dissipative unimodal maps) [14]. Albeit Ai vanishes, it is possible to express the divergence between trajectories with 
a form which conveniently generalizes the exponential function, namely the q-exponential form 

e (t) = exp^^ (Ag, t) (A,, > 0; < 1) , (2) 
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where exp^ x = [1 + {1 — q) x]^^^^ (exp^ x = e^), and Xq^ represents the generalised Lyapunov coefficient (the 
subscript s stands for sensitivity) [12]. Equation (2) can be looked as the solution of ^ = Xq^ S,'^' . Analogously to 
strong chaos (i.e., Ai > 0), if we optimize the entropy [15] Sq = ^ ^i^* i^i = Sbg) under the same type of 
constraint as before, we obtain p{x) oc expq{—(3qx), where (3q generalizes (3. 



II. THE CLASSICAL KICKED TOP MAP 



The classical kicked top corresponds to a map on the unit sphere x^ 
application 



1, corresponding to the following 



Xt+l = Zt 

yt+i = xt sin (a zt) + yt cos {a zt) 
zt+i = -X COS (a z) + yt sin (a zt) 



(3) 



where a denotes the kick strength. It is straighforward to verify that the determinant of the Jacobian matrix of (3) 

equals one, meaning that this map is conservative. It is therefore quite analogous to Hamiltonian conservative systems, 
the phase space of which consists of a mixing of regular (the famous KAM-tori [8]) and chaotic regions characterised, 
respectively, by a linear {qs = 0) and exponential {qs = 1) time evolution of ^ [12]. The region of separation presents 
a gs-exponcntial law for its sensitivity. In Fig. 1(a) we exhibit a trajectory where the various regions can be seen. 
In Figs. 1 (b,c,d) we see the time evolution of ^ for the three possible stages: regular, edge of chaos and chaotic, 
respectively. 
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FIG. 1. (a) Orbit of the a = 2.3 kicked top, where chaotic and regular regions are visible. The spherical phase space 
is projected onto x — z plane by multiplying the x and z coordinates of each point by Kjr where R — ^2 (1 — jjyj) and 
r = \/l — 3/^. (b-d) Time dependence of the sensitivity ^ to initial conditions (with ||Ar (0)|| = 10^'^'^) at (b) regular region 
(A; linear evolution), (c) edge of chaos (♦; qs -exponential evolution), and (d) chaotic region (T, ea;ponenttaZ evolution). 



It is worthy mentioning at this point that the quantum version of this map constitutes a paradigmatic example of 
quantum chaos. At its threshold to chaos, it has been verified a nonextensive behaviour (for details see [16]). 

We analysed here the sensitivity to initial conditions on the verge of chaos of (3), for several values of the kick 
strength a G [0, 4] averaged over a set of (typically 50) initial conditions for each value of a. More precisely, for fixed 
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a, aided by its typical orbits, wc determined a set of points in the regular-chaos border and then, for these points, 
determined the average value of ^ at time t. See typical results in Figs. 2 and 3. 
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FIG. 2. Time dependence of the averaged (over close to 50 points of the edge of chaos region) sensitivity (^) to initial 
conditions, for typical values of a. Insets: Same data but using a In^^- ordinate, where lnq(a;) = (a;^"'' — 1)/(1 — q) (Ini = In) . 
With this q-logarithm ordinate, the slope of the straight line is simply Aq, . 

We verify that the increase of a induces a gradual approach of to 1. This behaviour is in accordance with what 
was verified [12] for a non-linear system composed by two simplectically coupled standard maps. 
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FIG. 4. Representation of the orbits for typical values of a. As a increases we verify the emergence of chaotic regions that, 
for ac > — 3.2, fulfil the entire phase space. 

Summarising, we numerically analysed the sensitivity to initial conditions at the edge of chaos of the conservative 
classical kicked top, and found that its time evolution exhibits a gg-exponential behavior in all cases. For a = 0, the 
phase space is composed by a regular region, where the sensitivity depends linearly on time, hence qg = 0. As a 
increases, the top is more perturbed, hence chaotic regions emerge in phase space. Above some critical value Uc — 3.2, 
the chaotic region fulfils the entire phase space. Consistently, the usual exponential dependence (i.e., = 1) is 
recovered. These results can be useful to understand, within a nonextensive statistical mechanical framework, the 
everlasting metastable states that are known to exist in systems composed of many symplectically coupled maps [17], 
as well as in isolated many-body long-range- interacting classical Hamiltonians [18] . 

We would like to thank F. Baldovin and G.F.J. Ahahos for fruitful discussions, as well as FAPERJ, PRONEX and 
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